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SUMMARY 


The  work  reported  herein  demonstrates  the  feasibility 
of  applying  the  Euler  equations  to  calculating  flow  about 
supersonic  wing-body  combinations  at  high  angles  of  attack 
where  significant  effects  of  shock  waves  and  flow  separation 
occur.  Body  separation  vortices  and  leading-edge  separation 
from  the  wings  are  both  present.  It  is  necessary  to  input 
the  body  primary  separation  line  locations  to  apply  the  Euler 
code.  Special  data  correlations  to  yield  these  quantities 
were  carried  out. 

To  start  the  computations  in  a  supersonic  marching  code 
for  a  sharp-nosed  configuration,  it  was  necessary  to  have 
initial  solutions  for  circular  cones  with  separated  turbulent 
flow.  Such  solutions  were  developed  which  yield  conical 
pressure  distributions  in  good  agreement  with  data. 

The  code  was  applied  to  a  wing-body  combination  utilizing 
a  cone-cylinder  body  and  AR  *  1  delta  wings.  Calculations  were 
performed  for  Mach  numbers  up  to  5  and  angles  of  attack  up 
to  25°.  The  method  is  limited  in  angle  of  attack  only  by  the 
angle  of  shock  detachment  for  the  cone.  This  limit  could 
also  be  raised  if  starting  solutions  above  shock  detachment 
were  available.  The  body  pressure  distributions  are  generally 
in  good  agreement  with  data.  The  fin  normal  force  also  agrees 
with  experiment. 

A  special  comparison  was  made  between  the  load  induced 
on  the  fins  as  determined  by  the  Euler  code  and  as  calculated 
by  a  panel  method.  While  gross  loads  were  nearly  equal,  the 
spanwise  variation  of  the  pressure  differed. 


PRE'FACi 


This  technical  report  covers  the  work  performed  under 
Contract  N00014- 78-C-0490  from  December  1,  1980  to  November  30, 
1981.  The  program  is  sponsored  by  the  Office  of  Naval  Research, 
Arlington,  Virginia.  Dr.  Robert  E.  Whitehead  is  the  Scientific 
Officer, 
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1.  INTRODUCTION 


Wing-body  interference  at  high  angles  of  attack  and  high 
supersonic  speeds  is  characterized  by  significant  loss  of  fin 
lift  (Ref.  1)  and  by  the  lack  of  good  prediction  methods  to 
aid  missile  designers.  The  increased  emphasis  on  airbreathing 
missiles,  which  usually  utilize  noncircular  bodies,  has  further 
increased  the  need  for  a  reliable  and  rapid  predictive  method. 

For  the  low  angle  range  and  supersonic  Mach  numbers 
approaching  3,  panel  methods  (Ref.  2)  have  proved  useful.  At 
angles  above  the  linear  range  and  Mach  numbers  above  3,  the 
nonlinearities  due  to  compressibility  and  vorticity  require 
more  powerful  methods  for  their  predictions.  The  Navier- Stokes 
equations  represent  a  possibility  in  this  connection,  but  their 
use  for  three-dimensional  configurations  requires  large  amounts 
of  computer  time.  A  viable  alternative  appears  to  be  the 
Euler  equations  when  used  in  the  supersonic  marching  mode. 

This  mode  covers  much  of  the  range  of  angle  of  attack  and  Mach 
number  of  interest. 

The  Euler  equations  have  been  applied  to  supersonic 
missiles  by  investigators  at  the  Naval  Surface  Weapons  Center 
(Ref.  3)  and  by  the  present  investigators  (Refs.  1  and  4).  The 
former  investigators  treated  cases  of  attached  flow  while  the 
present  investigators  have  concentrated  on  separated  flows. 

In  Reference  1,  the  Euler  equations  were  applied  to  a  cruciform 
wing-body  combinations  at  zero  roll  angle  with  fins  having 
subsonic  leading  edges.  A  Kutta  condition  was  imposed  at  the 
sharp  leading  edges  to  make  the  computer  solution  unique, 
and  the  entire  flow  field  without  body  vortices  was  calculated. 
In  Reference  4,  a  special  study  was  made  of  generating  body 
vortices  with  the  Euler  code,  and  a  method  developed  which 
calculates  the  flow  field  well  except  near  secondary  separation. 
(The  method  can  include  secondary  separation  within  its  frame¬ 
work.)  The  primary  purpose  of  the  present  study  is  to  adapt 
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the  methods  previously  developed  to  the  calculation  of  wing- 
body  interference  with  the  effects  of  compressibility  and 
vortices  both  present.  Subsonic  leading  edges  are  treated. 

This  report,  after  discussing  the  general  background  of 
the  problem,  gives  a  r6sum6  of  the  new  calculation  method 
leaving  the  mathematical  details  to  the  earlier  reports.  Since 
the  location  of  the  primary  separation  line  is  an  input  to  the 
method,  the  report  next  addresses  this  question.  A  number  of 
calculative  examples  of  wing-body  interference  are  discussed 
to  illustrate  the  nature  of  the  phenomena  as  a  function  of 
angle  of  attack  and  Mach  number.  Finally,  the  Euler  method  is 
used  to  investigate  the  accuracy  of  panel  methods  for  calculat¬ 
ing  the  effect  of  body  vortices  on  fin  loading. 

The  work  herein  was  carried  out  under  Office  of  Naval 
Research  Contract  No.  N00014-78-C-0490  under  the  scientific 
direction  of  Dr.  Robert  E.  Whitehead. 


2.  GENERAL  CONSIDERATIONS 

The  present  investigation  is  not  meant  to  develop  a  pro¬ 
duction  computer  code  for  supersonic  cruciform  wing-body 
combinations,  but  is  meant  to  develop  the  necessary  methodology 
leading  to  such  a  code.  The  present  calculations  are  limited 
to  cruciform  wing-body  combinations  at  zero  roll  angle  although 
the  mathematical  steps  to  handle  non- zero  roll  angle  are  clear. 
The  Mach  number  range  is  supersonic  and  is  specifically  limited 
to  that  angle-of-attack  range  for  which  the  axial  Mach  number 
is  supersonic.  For  subsonic  axial  Mach  numbers,  problems  can 
be  solved  by  using  the  unsteady  Euler  equations  and  calculating 
hyperbolically  in  time  at  considerable  expense  in  computer  time 

At  the  present  time  sharp-nosed  bodies  of  revolution  are 
being  considered.  As  long  as  a  cone  solution  is  available  for 
the  given  angle  of  attack  and  Mach  number  to  start  the  solution 
it  should  proceed  subject  to  the  limitation  of  axial  supersonic 


6 


► 


Mach  number.  If  a  blunt  nose  is  used,  the  solution  can  proceed 
behind  any  station  of  the  starting  solution  for  which  the  axial 
Mach  number  exceeds  unity  everywhere.  For  noncircular  bodies 
of  revolution,  the  solution  can  still  be  carried  out  if  an 
appropriate  starting  solution  is  available  or  can  be  generated. 

In  the  present  work  any  pointed  body  of  revolution  can  be 
handled;  cone-cylinders,  ogive- cyl inder ,  etc.  However,  the 
wings  generally  have  thickness  distributions  of  a  kind  which 
correspond  to  blended  wing-body  junctures.  This  results  from 
the  way  in  which  conformal  mapping  is  used  in  the  method. 

The  flow  around  wing-body  combinations  at  supersonic  speed 
and  moderate  to  high  angles  of  attack  is  dominated  by  compress¬ 
ibility  and  vortex  nonlinearities.  The  Navier- Stokes  equations 
are  applicable  to  these  flows  but  take  much  computer  time.  The 
application  of  the  Euler  equations  to  these  problems  is  moti¬ 
vated  by  the  possibility  of  getting  faster  solutions  and  by 
the  fact  that  diffusion  of  vorticity  not  accounted  for  by  the 
Euler  equations  is  of  little  importance  in  many  problems  of 
missile  aerodynamics.  The  convection  of  vorticity  as  in 
trailing  vortices  is  of  importance,  and  can  be  accounted  for 
by  the  Euler  equations  provided  the  vorticity  is  properly 
introduced  into  the  flow.  The  vorticity  generated  at  curved 
shock  waves  is  generated  external  to  the  solid  surfaces  and 
appears  in  the  Euler  solutions.  The  vorticity  which  arises  at 
the  surface  such  as  body  vortices  or  leading-edge  vortices,  must 
be  controlled  in  the  Euler  equations  by  use  of  appropriate 
boundary  conditions. 

The  conditions  required  to  control  the  vort'city  discharged 
into  the  flow  at  solid  boundaries  include  a  knowledge  of  the 
separation  line  locations  together  with  additional  conditions 
required  to  tie  down  all  the  dependent  variables.  In  the  first 
year's  work,  we  showed  how  the  Kutta  condition  at  a  subsonic 
leading  edge  plus  certain  auxiliary  boundary  conditions  were 
sufficient  to  allow  calculation  of  separated  flows  about  wings 
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alone  and  the  fins  of  wing-body  combinations.  The  problem 
here  is  simplified  by  the  fact  that  the  separation  location  at 
the  sharp  leading  edge  is  known. 

For  bodies  the  locations  of  the  separation  lines  at  super¬ 
sonic  speeds  are  not  well  known.  We  will  subsequently  enlarge 
on  this  subject  as  an  aid  in  applying  the  Euler  equations  to 
the  problem.  However,  the  position  of  the  separation  line  is 
not  enough  information  to  tie  down  all  the  dependent  variables. 
Certain  other  conditions  relating  to  the  direction  of  the 
vortex  lines  leaving  the  separation  line  were  used  together  with 
mean  value  conditions  for  the  pressure  and  density  in  the  separa 
tion  sheet  at  separation.  These  conditions  were  sufficient  to 
determine  the  rate  of  vorticity  shedding  at  the  separation  point 
as  shown  by  comparison  (Fig.  22,  Ref.  4)  with  the  data  of 
Oberkampf  (Ref.  5).  A  parametric  study  of  the  boundary  condi¬ 
tions  revealed  a  general  insensitivity  of  the  flow  field  to 
the  directions  of  the  vortex  sheets  leaving  the  separation  line. 


3.  RESUME  OF  CALCULATIVE  METHOD 

A  complete  description  of  the  mathematical  method  is 
lengthy  and  is  contained  in  References  1  and  4.  However,  a 
rdsumd  of  the  method  is  given  here  as  an  aid  to  the  future 
sections  of  the  report. 

3.1  Coordinates  and  Mesh  Generation 

The  physical  coordinate  system  is  an  x,y,z  Cartesian  system 
as  shown  in  Figure  1  for  a  cone-cylinder  combination.  The 
origin  is  at  the  apex  of  the  cone  with  the  z-coordinate  rear¬ 
ward.  The  y  coordinate  is  vertically  upward.  At  each  section 
the  configuration  has  a  cross  section  which  is  mapped  con¬ 
formally  into  the  unit  circle  in  the  transformed  space.  This 
transforms  the  body  or  wing  body  into  a  cylinder.  The  bow 
wave  is  generally  not  concentric  with  the  body.  Accordingly, 
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with  the  cylinder  fixed,  a  radial  linear  transformation  which 
depends  on  <J>  is  applied  to  make  the  bow  wave  concentric  with 
the  body. 
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(1) 


This  makes  £;  =  0  on  the  body  and  unity  on  the  bow  wave.  An 
orthogonal  mesh  can  now  be  laid  out  and  carried  back  to  the 
transformed  space  and  physical  space  where  it  is  no  longer 
orthogonal. 

Under  these  transformations,  the  metric  terms  are  obtained 
algebraically  from  the  transformation  equations.  For  a  wing- 
body  cross  section  the  Joukowski  transformation  carries  it 
into  a  circle  as  shown.  The  Jacobian  of  the  transformation 
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shows  singularities  at  the  tip  of  a  sharp  wing  or  at  the  wing- 
body  juncture.  To  avoid  the  singularity  at  the  wing-body 
juncture,  the  contour  of  the  configuration  has  been  taken  as 
a  circle  of  radius  1  +  e,  where  e  is  a  small  number.  The 
resulting  wing-body  combination  is  blended  as  shown  in  the 
following  sketch; 


The  singularity  at  the  leading  edge  of  the  fin  is  handled 
by  a  Kutta  condition  as  subsequently  described. 


3.2  Governing  Equations 

The  governing  equations  for  the  steady  inviscid  flow  of 
an  ideal  compressible  fluid  in  three  dimensions,  are  the  Euler 
equations.  In  Cartesian  coordinates  the  equations  are 


it  +  il  + 

ax  3y  3z 


0 


A 

where  E, 


F,  and  G  are  defined  as 


(2) 


pu 

pv 

pw 

2 

pu  +  kp 

puv 

puw 

F  = 

G  = 

puv 

2 

pv  +  kp 

pvw 

puw 

pvw 

2 

pw  +  pk 

The  equations  involve  five  dependent  variables,  three  velocity 
components,  u,  v,  and  w  and  two  thermodynamic  variables  of 
state,  p  andp.  Another  relationship  is  provided  by  the  con¬ 
stancy  of  total  temperature  or  total  enthalpy.  These  equations 
in  the  transformed  space  are  still  fully  conservative,  but 
the  radial  stretching  operation  to  make  the  bow  shock  correspond 
to  £  =  1  makes  the  equations  quasi-conservative. 

3.3  Boundary  Conditions  at  Bow  Wave  and  Regular 
Points  on  the  Solid  Boundary 

At  the  bow  wave,  the  Rankine- Hugoniot  relationship  must  be 
satisfied.  By  utilizing  information  for  data  at  the  shock  and 
one  point  below  the  shock,  a  prediction-correction  scheme  is 
used  to  determine  the  shock  shape  as  the  calculation  proceeds 
axially. 

The  flow  normal  to  the  wing-body  surface  at  all  nonsingular 
points  must  be  zero  for  steady  flow.  The  values  of  the 
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dependent  variables  are  calculated  at  the  next  station  on  the 
body  using  the  MacCormack  prediction-correction  scheme  used  for 
the  body  tangency  condition.  A  scheme  due  to  Abbett  (Ref.  6) 
is  used  to  correct  the  boundary  values.  In  this  prediction- 
correction  scheme,  simple  expansion  or  compression  waves  are  used 
at  the  surface  in  the  tangent  plane.  The  possibility  of  cor¬ 
recting  the  surface  values  of  the  entropy  for  shock  waves  which 
reach  to  the  body  surface  is  included.  The  boundary  conditions 
are  applied  at  the  wing  and  body  surfaces. 

3.4  Boundary  Conditions  at  a  Sharp  Subsonic  Leading  Edge 

At  a  wing  tip  the  metric  can  be  expanded  in  powers  of 
small  perturbations  of  the  radial  and  angular  variables  from 
their  values  at  the  wing  tip.  Using  those  results  to  expand 
the  Euler  equations  in  the  neighborhood  of  the  wing  tip,  it 
is  found  (Ref.  1)  that  the  equations  are  indeterminate. 
Accordingly,  the  Euler  equations  are  satisfied  by  any  con¬ 
tinuous  values  of  the  variables  approaching  the  wing  tip.  It 
was  also  found  that  the  same  is  true  of  jump  conditions  across 
any  sheet  of  discontinuity  at  the  tip.  Accordingly,  we  have 
five  independent  conditions  that  can  be  specified  at  the  sub¬ 
sonic  edge.  The  three  conditions  which  we  have  chosen  to  use  on 
the  basis  of  physical  arguments  are: 

1.  The  static  pressure  is  continuous  across  the  sheet. 

2.  The  vortex  sheet  leaves  the  leading  edge  in  the  y  =  0 
plane  (Kutta  condition). 

3.  The  total  enthalpy  is  constant. 

The  other  two  conditions  which  have  been  specified  are: 

4.  The  density  at  the  wing  tip  is  the  average  of  the 
density  on  the  wing  on  each  side  of  the  sheet. 

5.  The  initial  direction  of  the  vortex  lines  in  the 
sheet  is  the  average  of  the  boundary  streamlines 
just  inboard  of  the  tip  on  the  upper  and  lower 
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surfaces  provided  it  is  outboard.  Otherwise  it  is 
taken  tangent  to  the  wing  leading  edge. 

These  conditions  are  adequate  to  remove  the  indetermin¬ 
ateness  of  the  Euler  equations  at  the  leading  edge.  The 
significance  of  each  of  these  conditions  in  determining  the 
gross  flow  is  another  matter. 

3.S  Boundary  Conditions  at  Separation  Lines 
on  a  Circular  Body 

At  a  separation  line  on  a  circular  or  noncircular  body, 
the  Euler  equations  are  not  indeterminate  as  at  a  sharp 
subsonic  leading  edge.  However,  the  tangency  condition  is 
not  valid.  The  flow  variables,  u,  v,  w,  p,  and  p  must  satisfy 
the  conditions  for  a  tangential  discontinuity.  What  we  have 
done  is  let  the  Euler  equations  capture  the  discontinuities 
which  are  determined  by  global  conditions,  not  local  condi¬ 
tions.  We  rather  express  the  mean  values  in  the  sheet  based 
on  various  physical  assumptions. 

1.  The  pressure  across  the  sheet  is  continuous.  This 

is  expressed  by  the  equation  that  the  static  pressure 
at  separation  is  the  average  of  the  body  pressure 
above  and  below  the  separation  line. 

2.  The  mean  density  within  the  sheet  is  specified  by  a 
relationship  similar  to  that  for  the  pressure. 

3.  The  total  enthalpy  is  constant. 

In  addition  to  these  conditions,  certain  conditions  on 
the  position  and  direction  of  the  vortex  lines  at  separation 
are  also  needed  to  close  the  problem.  Figure  2  shows  three 
angles  which  specify  this  information  for  the  vortex  lines 
leaving  the  surface  of  a  round  body.  The  angle  <j>s  specifies 
the  polar  angle  of  the  separation  line  in  the  crossflow  plane. 
The  angle  <j>c  is  the  angle  of  the  plane  containing  the  separa¬ 
tion  line  and  the  velocity  vector  with  respect  to  the  tangent 
plane.  The  angle  <|>a  is  related  to  the  sweepback  of  the 
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velocity  vector,  being  the  complement  of  the  sweepback  angle. 

The  velocities  uQ,  vQ,  and  wQ  at  the  separation  point  are: 

uo  =  qosin^acos(0s  -  *c) 

vQ  =  q0sin4»asin(^s  -  ♦c)  (3) 

wo  =  q0cos^a 

In  these  equations,  the  two  angles  <|>  and  (4>  -  <p  )  yield  the 

a  S  L 

ratios  u  /q  ,  v  /q  ,  and  w  /q  .  The  three  conditions  above 
yield  qQ  so  that  uQ,  vQ,  and  wq  are  known  from  Equation  (3). 

The  first  two  conditions  yield  pQ  and  pQ.  Thus,  five  condi¬ 
tions  yield  the  five  dependent  variables  at  separation. 

A  parametric  study  (Ref.  4)  of  the  flow  dependence  on  <J>  , 

a 

<t>c,  and  <t>s  indicates  that  the  rate  of  vortex  shedding  and  the 
height  of  the  separation  region  are  insensitive  to  the  values  of 
the  first  two  parameters.  The  important  parameter  is  4>s. 

3.6  Initial  Conditions 

It  is  necessary  to  initialize  the  program  at  some  axial 
station  at  which  a  crossflow  solution  is  known.  The  common 
case  is  that  of  a  pointed  body  of  revolution  which  is  initial¬ 
ized  by  the  solution  for  the  tangent  cone  to  the  nose  at  the 
given  Mach  number  and  angle  of  attack.  If  a  cone  solution 
close  to  the  desired  solution  is  known,  it  can  be  used  to 
determine  the  precise  starting  solution.  If  a  certain  cone 
half  angle  is  needed,  for  instance,  we  can  start  the  solution 
with  a  known  half  angle  and  then  change  to  a  different  half 
angle  as  soon  as  the  calculation  is  started.  The  solution  is 
then  continued  downstream  until  the  effect  of  the  initial 
solution  has  died  out.  This  is  determined  by  testing  the 
solution  for  conical  flow  as  it  proceeds  downstream.  This 
procedure  is  valid  up  to  the  shock-detachment  angle  of  attack. 
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The  present  program  can  also  be  used  to  study  wing-alone 
characteristics  at  high  angles  of  attack.  The  solution  is 
started  with  a  known  solution  for  a  cone,  but  the  cross  section 
is  changed  from  a  circle  to  an  ellipse  of  high  eccentricity  as 
the  calculation  proceeds  downstream  to  form  a  wing.  Again  the 
solution  is  carried  downstream  until  it  is  conical  to  generate 
a  starting  solution  for  the  real  wing.  While  not  all  wings 
of  interest  can  be  handled  this  way,  many  can. 

In  some  cases  bodies  have  blunt  noses  which  cause  detached 
bow  waves  and  imbedded  subsonic  flow  regions.  In  this  case  the 
starting  solution  can  be  generated  by  an  unsteady  Euler  code 
(Ref.  7) ,  and  the  marching  program  can  be  engaged  downstream  of 
the  region  of  imbedded  subsonic  flow.  A  question  arises  con¬ 
cerning  how  to  start  a  solution  for  a  body  with  a  conical  nose 
when  separation  starts  at  the  apex.  This  point  will  be  sub¬ 
sequently  addressed. 

4.  PRIMARY  SEPARATION  LINE  LOCATIONS 
4.1  Circular  Cylinders 

The  location  of  the  separation  point  for  circular 
cylinders  is  of  interest  for  what  light  it  might  shed  on  the 
location  of  the  separation  point  on  the  cylindrical  part  of 
a  body  of  revolution.  One  does  not  have  to  go  very  deeply 
into  the  subject  to  appreciate  the  complexity  of  the  problem 
for  the  circular  cylinder.  In  this  connection  consider 
Figure  3  from  Reference  8  which  shows  the  separation  positions 
on  smooth  circular  cylinders  as  a  function  of  Reynolds  number 
for  incompressible  flow.  In  the  subcritical  region  the  boundary 
layer  is  laminar  at  the  separation  point  with  no  reattachment. 

In  the  critical  region  a  laminar  separation  bubble  forms  with 
turbulent  reattachment  and  subsequent  turbulent  separation. 

In  the  supercritical  region  the  transition  occurs  in  front  of 
separation  and  the  separation  is  all  turbulent.  The  separation 
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locations  were  determined  from  local  skin- friction  measure¬ 
ments  . 

The  effect  of  Mach  number  on  the  location  of  the  separa¬ 
tion  point  is  of  interest  as  it  might  apply  to  estimation  of 
primary  separation  lines  on  bodies  of  revolution  at  angles  of 
attack.  The  increase  of  Mach  number  up  to  about  M^  =  0.4  will 
not  cause  any  part  of  the  flow  to  exceed  the  speed  of  sound, 
but  above  this  Mach  number  inviscid  flow  theory  predicts 
imbedded  regions  of  supersonic  flow  and  shock  waves.  It 
could  be  expected  that  the  shock  waves  would  have  an  effect 
on  separation,  and  if  strong,  may  control  the  position  of 
separation. 

Very  little  data  are  available  on  this  point,  the  recent 
work  (1977)  of  Murthy  and  Rose  (Ref.  9)  being  a  notable  excep¬ 
tion.  Some  interesting  facts  can  be  gleaned  from  a  study  of 
these  data,  which  include  the  critical  and  supercritical 
regimes.  The  ranges  of  the  test  paramters  are: 

.083  x  105  <  ReD  <  0.5  x  106 

0.25  <  M  <1.2 
—  00  — 

The  Reynolds  number  is  based  on  cylinder  diameter.  For 
M^  >  0 . 5  and  ReD  =  0.5  x  106  boundary- layer  transition  occurs 
on  the  front  of  the  cylinder,  and  the  separation  appears  to 
be  turbulent.  It  is  of  interest  that  this  supercritical  state 
is  achieved  at  a  Reynolds  number  for  which  the  data  of 
Achenbach  (Fig.  3)  are  in  the  critical  regime.  This  result 
may  not  be  entirely  due  to  compressibility. 

The  following  separation  points  were  obtained  (Ref.  9) 
from  the  skin-friction  measurements  made  by  a  buried  wire  with 
supercritical  crossflow. 
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These  points,  shown  in  Figure  4  as  a  function  of  Mach  number, 
indicate  a  strong  effect  of  Mach  number.  In  Figure  5  a  plot 
is  shown  of  the  skin  friction  and  pressure  distribution  around 
the  cylinder  for  the  M^  =  0.8  case.  The  data  exhibit  the 
interesting  phenomenon  that  the  pressure  becomes  uniform  for 
ip  >  75°,  but  separation  is  not  seen  until  <J>  ~  80°.  An  angular 
separation  between  minimum  pressure  onset  and  separation  of 
about  6°  -  7°  occurs  at  nearly  constant  pressure.  There  is  no 
pressure  rise  in  the  experimental  pressure  distribution  before 
separation  which  can  be  used  to  estimate  the  separation  point 
by  such  criteria  as  that  due  to  Stratford  (Ref.  10).  It 
appears  that  the  separation  position  may  be  controlled  by  some 
other  mechanism  than  an  adverse  pressure  gradient  in  the 
Stratford  sense. 

Another  noteworthy  fact  is  that  there  is  no  evidence  of 
a  shock  wave  in  any  of  the  surface  pressure  distributions 
although  a  weak  one  might  exist.  Somehow  the  boundary- layer 
separation  must  act  to  eliminate  the  shock  pressure  discon¬ 
tinuity  at  the  surface.  This  raises  the  conjecture  of  whether 
separation  occurs  at  or  near  the  shock  location.  Let  us 
examine  this  question  by  comparing  theoretical  and  experimental 
results.  The  use  of  potential  theory  and  CFD  calculations 
will  yield  a  position  of  the  theoretical  shock  position.  Such 
a  calculation  for  M  =  0.51  is  given  in  Reference  11.  This 
point  is  shown  in  Figure  4  for  comparison  with  the  experimental 
separation-point  locations.  It  is  clearly  well  behind  the 


probable  position  of  the  separation  point.  This  is  not 
surprising  in  view  of  the  fact  that  the  full  potential  equa¬ 
tion  is  known  to  calculate  shock  locations  behind  the 
experimental  ones  for  airfoils  in  transonic  flow. 

We  will  address  the  question  of  the  application  of  the 
circular  cylinder  results  to  the  cylindrical  portion  of 
todies  of  revolution  at  angle  of  attack  subsequently. 

4.2  Pointed  Cones 

4.2.1  Preliminary  remarks.  If  a  cone  is  pitched  up  to 
an  angle  of  attack  equal  to  its  half  angle,  0  ,  the  top 
meridian  is  parallel  to  the  free  stream  and  the  local  pressure 
is  close  to  the  free-stream  pressure.  Further  increases  in 
angle  of  attack  are  usually  accompanied  by  separation  on  the 
leeward  side  of  the  cone  starting  at  the  top  meridian  and 
moving  downward  as  the  angle  of  attack  is  increased.  It  is 
observed  experimentally  that  flow  separation  occurs  along  a 
straight  line  at  a  constant  value  of  azimuth  angle,  4>s .  The 
flow  pattern  on  the  attached  side  of  the  cone  appears  conical 
even  for  a  finite  cone.  It  is  known  that  the  position  4>s 
depends  on  the  ratio  a/0c  to  a  significant  degree. 

Another  parameter  that  might  be  expected  to  have  a  signi¬ 
ficant  effect  on  the  value  of  <f>s  is  the  crossflow  Mach  number. 

If  we  draw  an  analogy  between  the  flow  in  the  crossflow  plane 
of  a  cone  and  that  for  a  circular  cylinder  at  the  same  cross- 
flow  Mach  number,  we  might  expect  a  change  in  the  qualitative 
nature  of  the  flow  with  increases  in  Mach  number  when  super¬ 
sonic  flow  occurs  in  the  crossflow  plane.  For  this  reason  the 
data  for  Mn  greater  than  a  certain  value  might  be  expected  to 
correlate  differently  than  that  for  Mach  numbers  less  than  this 
value. 

An  interesting  discussion  of  the  flow  over  cones,  including 
separated  flow,  treating  both  experiment  and  theory  is  given 
by  Peak  and  Toback  in  Reference  12.  In  that  reference  the 
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authors  have  collected  together  some  of  the  data  on  primary 

cone  separation  angle,  $  .  Both  laminar  and  turbulent  flow 

s  i 

are  represented  as  well  as  low  and  high  crossflow  Mach  numbers. 
We  will  now  discuss  the  available  data  before  presenting  the 
correlations . 

4.2.2  Description  of  available  data.  We  will  first 
consider  the  data  on  primary  separation  for  cones  from  a  number 
of  references.  Probably  the  best  collection  of  data  is  that 
of  Rainbird  in  Reference  13,  Figure  14.  Rainbird  tested  cones 
with  0  =  5°  at  M  =  1.8  and  4.25  over  a  wide  angle- of - attack 

C  00 

range.  The  measurements  were  made  in  the  NAE  5- Ft  Trisonic 
Blowdown  Tunnel  on  a  54- inch  long  cone.  The  turbulence  and 
noise  of  the  tunnel  induced  early  transition  within  10  percent 
of  the  body  length  from  the  nose.  The  cone  had  essentially 
the  same  temperature  as  the  tunnel  total  temperature  so  that 
little  heat  transfer  occurred.  The  location  of  the  separation 
and  attachment  lines  was  determined  using  the  oil-film 
technique.  Results  on  a  40.6-inch  long  cone  of  0c  =  12.5°  are 
also  available  from  tests  under  similar  conditions.  In 
Reference  14  the  data  for  the  12.5°  cone  are  presented,  but  no 
further  primary  separation  data  are  included  not  already  con¬ 
tained  in  Reference  13.  The  results  are  given  in  Table  1. 

Another  useful  set  of  data  is  given  by  Nebbeling  and 
Bannink  in  Reference  15.  Their  cone  of  half  angle  7.5°  was 
tested  in  two  supersonic  blowdown  wind  tunnels  at  the  University 
of  Technology,  Delft.  The  0.27m  by  0.27m  tunnel  was  used  to 
provide  flow-surface  visualization.  A  mixture  of  oil,  detergent 
and  titanium  dioxide  was  used  as  a  surface  film  to  provide  the 
flow  visualization.  The  tests  were  conducted  at  =  2.94  at  a 
Reynolds  number  per  meter  of  5  *  10  .  Data  on  the  secondary 
separation  line  position,  4>S2>  are  also  given.  The  results 
agree  with  the  turbulent  flow  results  of  Rainbird  (Ref.  13). 

The  model  was  cooled  down  during  the  runs  but  it  is  thought  that 
the  heat  transfer  rate  was  small.  The  flow  visualization 
pictures  are  of  excellent  detail  and  quality. 
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McElderry  (Ref.  16)  provides  some  data  on  separation  line 
position  at  =  6.0.  His  tests  were  performed  in  the  50-inch 
continuous  flow  hypersonic  wind  tunnel  of  the  AEDC  von  KdrmSn 
facility  with  pointed  and  blunted  cones.  The  nominal  test 
conditions  were:  =  6.05,  Re/ft  =  5  x  10^,  PQ  =  280  psia, 

Tq  =  850°R.  The  model  surface  temperature  started  out  well 
below  the  tunnel  stagnation  temperature.  Based  on  a  model 
length  of  38.06  in.,  the  Reynolds  number  based  on  length  is 
about  1.27x10  .  The  flow  visualization  was  carried  out  using 
a  white  pigmented  silicone  oil  on  the  dark  model.  The  model  was 
inserted  into  the  airstream  for  a  few  minutes  and  then  photo¬ 
graphed  from  top  and  side  viewing  ports.  Data  are  given  for 
the  primary  and  secondary  line  location  for  a  blunted  cone. 

The  data  are  for  the  turbulent  separation  region  which  occurred 
behind  the  tip  of  the  cone.  There  is  a  band  of  <f>  given  in 
Figure  12  of  Reference  16  which  is  included  in  the  results 
of  Table  1.  It  was  noted  that  at  a/0c  of  unity,  flow  separa¬ 
tion  had  not  appeared. 

Peake,  Fisher,  and  McRae  (Ref.  17)  have  conducted  a 
combined  wind  tunnel,  flight,  and  analytical  investigation  of 
the  flow  about  a  cone  of  5°  semiapex  angle.  The  wind-tunnel 
tests  were  conducted  in  the  Ames  Research  Center  6-  by  6-Foot 
Supersonic  Wind  Tunnel.  The  flight  tests  were  conducted  with 
the  cone  mounted  on  a  nose  boom  on  an  F-15  airplane.  The 
cone  was  44.5-inches  long  with  roughness  at  4.9  inches  from 
the  apex  to  ensure  early  transition.  Since  oil  film  techniques 
were  not  appropriate  to  determine  separation  line  locations, 
obstacle  blocks  were  used.  The  obstacle  blocks  were  similar 
to  a  total-head  tube  flush  with  the  cone  surface.  The  measured 
results  are  given  in  Table  1.  The  correspondence  between  the 
flight  results  and  the  wind-tunnel  results  is  remarkable. 

Feldhuhn  and  Pasiuk  (Ref.  18)  have  also  provided  data  on 
primary  and  secondary  separation- line  position  for  a  cone  of 
5°  semiapex  angle  at  a  Mach  number  of  about  6.0.  Their  tests 
were  conducted  in  the  hypersonic  tunnel  of  the  old  Naval 
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Ordnance  Laboratory.  The  cone  was  sting  supported  and  had  a 
slightly  blunted  nose  of  radius  .0313  in.  compared  to  a  base 
diameter  of  1.25  in.  The  overall  length  was  13.967  in.  The 
separation- line  locations  were  determined  using  flow  visualiza¬ 
tion  based  on  a  mixture  of  lamp  black  and  silicone  oil.  After 
application  of  the  mixture,  the  model  was  shielded  while  the 
tunnel  was  started.  The  model  was  exposed  to  the  flow  for 
about  30  seconds  and  then  was  again  shielded.  Measurements 
of  the  angular  position  of  the  separation  lines  were  made 
outside  the  tunnel  at  length  stations  of  5.7,  8.5,  and  11.4  in. 
The  angular  positions  are  accurate  to  about  ±5°.  Secondary 
separation  was  found  for  angles  of  attack  of  20°  or  less. 

The  data  are  given  in  Figure  25  of  Reference  18.  The  data  are 
for  a  highly- cooled  cone. 

4.2.3  Correlation  of  the  data.  A  correlation  of  the 
data  of  Table  1  was  developed  from  the  observation  that  the 
static  pressures  on  cones  for  the  supersonic  case  can  be  shown 
to  follow  an  hypersonic  similarity  law  of  the  following  type. 


P  = 


1/2  pV^sin^a 


(4) 


Neglecting  differences  between  the  sine  of  an  angle  and 
the  angle  itself.  Equation  (4)  can  be  changed  to  an  equivalent 
form 
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where  F  is  a  new  function.  The  quantity  qmsin  a  is  the  com¬ 
ponent  of  qOT  in  a  direction  normal  to  the  cone  axis  within  the 
small  angle  approximation.  If  a  cone  of  small  semiapex  angle 
at  angle  of  attack  can  be  thought  of  as  a  swept  wire,  then  the 
simple-sweep  laminar  flow  analogy  of  R.T.  Jones,  Reference  19, 
can  be  invoked.  In  this  case  the  flow  in  the  crossflow  plane 
completely  determines  the  viscous  flow.  Assuming  that  the 


form  of  the  azimuthal  pressure  distribution  still  follows 
Equation  (4)  even  with  separated  flow,  then  the  primary  separa¬ 
tion  point  location,  <j>  ,  should  depend  only  on  a/8c  and  Mn. 

This  functional  relationship  thus  suggests  the  form  of  a  trial 
correlation.  Another  noteworthy  point  is  that  generally 
separation  on  a  cone  starts  approximately  at  a  =  9C.  If  <J>S1 
is  plotted  versus  a/0c,  the  lines  for  different  values  of  Mn 
should  thus  have  a  common  point  near  a/0c  =  1 

With  regard  to  correlating  the  separation  angle  4>s  versus 
a/0c  for  various  certain  statements  can  be  made.  It  is  sus¬ 
pected  that  for  subcritical  crossflow  Mach  numbers  the  data  should 
correlate  differently  from  those  for  supercritical  crossflow 
Mach  numbers.  This  turns  out  to  be  the  case  if  Mn  <  0.6.  For 
the  transonic  range  taken  to  be  0.6  <  Mn  <  1.2  deviations  from 
the  correlation  for  Mn  <  0.6  would  be  anticipated.  The 
0.6  <  Mn  <  1.2  correlation  might  be  expected  to  hold  for 
Mn  >  1.2  since  the  separation  point  is  now  controlled  by  cross- 
flow  shocks  rather  than  by  smooth  pressure  rise  as  for  Mn  <  0.6. 
Accordingly,  we  now  have  correlated  the  data  differentiating 
three  ranges  of  Mach  number 
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.2  < 

Mn 

The  resultant  correlation  is  shown  in  Figure  6.  The  symbols  are 
defined  in  Table  1. 

Figure  6  exhibits  a  number  of  interesting  points.  Con¬ 
sider  the  circles  for  a  12.5°  cone  and  the  diamonds  for  a 
7.5  degree  cone.  They  both  exhibit  a  sharp  decrease  in  <f>s^  as 

M  becomes  greater  than  M„  =  0.6.  All  the  data  correlate  on 
n  n 

two  lines  as  shown  within  fairly  good  accuracy.  (This  does  not 
prove  that  all  other  data  will  fall  into  the  correlation  since 
tests  varying  Mn  and  a/0c  systematically  are  not  available.) 

As  Mn  exceeds  0.6  the  value  of  s ^  falls  from  the  upper  curve 
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to  the  lower  curve.  While  the  correlation  shows  a  discon¬ 
tinuous  change  at  Mn  =  0.6,  the  actual  change  is  steep  but 
continuous . 

It  is  noted  that  <J>S^  does  not  approach  180°  at  a/0c  =  1 
but  a  small  region  of  separated  flow  is  indicated.  Questions 
will  arise  for  cones  of  large  0c  when  M  exceeds  0.6  for  values 
of  a/0c  less  than  unity.  Our  data  do  not  span  this  regime. 

We  now  consider  how  the  foregoing  information  can  be  used 
in  predicting  the  primary  separation  lines  for  cone-cylinder 
combinations . 

4.2.4  Applicable  cone-cylinder  results.  In  addition  to 
the  data  already  discussed,  a  systematic  series  of  oil-flow 
pictures  is  available  for  a  cone-cylinder  combination  for  a 
wide  range  of  angles  of  attack  and  Mach  numbers.  The  investi¬ 
gation  was  carried  out  by  Landrum  and  Babb  (Ref.  20)  on  a 
9.47°  half-angle  cone  followed  by  a  cylinder  3-2/3  diameters 
long.  The  angle  of  attack  range  was  -4°  to  60°  and  the  Mach 
number  range  was  1.6  to  4.63.  Oil-flow  pictures  are  given 
for  the  entire  range  of  angle  of  attack  and  Mach  number  from 
which  values  of  the  cone  primary  separation  angle  can  be 
measured.  An  index  line  for  <J>  =  90°  helped  greatly  in  this 
process,  but  nevertheless  the  accuracy  of  the  estimated 
angles  of  d> s ^ ,  are  estimated  to  be  ±7.5°.  The  method  of 
estimating  the  angle  of  the  primary  separation  line  from  the 
photographs  of  the  oil  flow  pattern  is  described  in  Appendix  A. 
The  positions  so  obtained  are  summarized  in  Table  2. 

Table  2  gives  the  measured  separation  line  positions 
both  on  the  conical  nose  and  on  the  cylindrical  part  of  the 
cone-cylinder  body  over  wide  Mach  number  and  angle-of-attack 
ranges.  We  now  consider  the  measured  angle  near  the  base  of 
the  body  at  x/l  =  0.925  and  compare  these  with  the  angles  for 
circular  cylinders  in  normal  crossflow  in  Figure  4.  Following 
Figure  4  and  the  correlation  of  Figure  31  in  Reference  21,  we 
have  plotted  <ps  versus  crossflow  Mach  number  Mn  for  all  the 
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data,  and  the  results  are  presented  in  Figure  7.  As  Mn 
increases  from  a  subsonic  value,  the  value  of  increases 

and  becomes  fairly  constant  above  Mn  =  1.2  at  a  value  of 
about  110°.  The  shape  of  the  correlation  is  similar  to  that 
of  Figure  31,  Reference  21,  since  the  data  for  the  two  correla¬ 
tions  overlap.  However,  the  present  correlation  goes  up  to 
Mn  =  4.0  while  the  former  goes  up  only  to  Mn  =  1.9.  It  seems 
clear  that  the  data  from  the  circular  cylinder  in  normal 
crossflow  and  the  cone-cylinder  bodies  at  angle  of  attack 
correlate  well  together. 

As  a  matter  of  interest,  let  us  see  how  the  correlation 
curves  of  Figures  6  and  7  predict  the  primary  separation  angle 
locations  on  the  cone-cylinder  body  of  Reference  20  on  the 
cone  at  x/l  =  0.325  and  on  the  cylindrical  part  of  the  body  at 
x/l  =  0.925.  Since  the  data  from  these  runs  have  been  used 
in  the  correlation  of  Figure  7,  we  expect  accuracy  of  predic¬ 
tion  at  x/l  =  0.925  to  be  within  the  scatter  of  correlation. 
However,  the  prediction  of  x/l  =  0.325  is  based  solely  on 
cone-alone  data.  The  comparisons  are  made  in  Figure  8  for 
four  Mach  numbers.  Examining  first  the  comparisons  for 
x/l  =  0.925,  we  see  that  the  agreement  is  within  the  scatter 
of  Figure  7  as  expected.  Examining  now  the  comparisons  for 
x/l  -  0.325,  we  see  the  agreement  is  almost  as  good  as  for 
x/l  =  0.925  except  at  =  1.60  at  a  =  20°.  This  condition  is 
just  before  the  onset  of  the  transonic  regime  at  Mn  =  0.6  as 
shown  on  the  abscissa.  The  rapid  fall  in  $s  for  a  cone  alone 
at  this  Mach  number  could  explain  this  disagreement. 


5.  CALCULATED  RESULTS  FOR  POINTED  CONES 
5.1  Introductory  Remarks 

To  start  the  calculations  for  most  supersonic  missile 
configurations,  a  sharp  cone  is  fitted  to  the  nose  and  the  cone 
solution  is  used  for  the  initial  conditions.  For  sharp  cones 
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up  to  an  angle  of  attack  approximately  equal  to  the  cone  semi¬ 
apex  angle,  the  flow  is  attached  and  no  particular  problem 
exists  in  determining  the  starting  solution.  However,  for 
angles  of  attack  above  the  semiapex  angle,  the  flow  is 
usually  separated  and  a  question  arises  concerning  the  proper 
starting  solution.  We  have  observed  in  the  previous  section 
that  experimentally  separated  flows  on  cones  are  essentially 
conical.  What  we  now  consider  is  the  starting  solution  for 
the  cone  with  a/6c  greater  than  unity  with  and  without  a  Kutta 
condition  imposed. 

5.2  Starting  Solution  with  No  Kutta  Condition 

The  calculation  was  made  for  a  conical  nose  of  fineness 
ratio  3  so  that  the  semiapex  angle  is  9.465°.  The  free-stream 
Mach  number  is  3.0  and  the  angle  of  attack  is  25°.  The  first 
problem  is  to  start  the  solution  since  the  initial  conditions 
required  are  the  flow  field  in  a  crossflow  plane  at  some  small 
value  of  axial  distance,  z.  The  solution  is  continued  until 
the  effect  of  the  starting  solution  dies  out.  For  conical 
flow  we  can  test  the  conicity  of  the  flow  to  determine  when  the 
effect  of  incorrect  initial  conditions  are  damped  out. 

Since  there  is  no  length  scale  in  the  problem,  arbitrary 
units  of  axial  length  have  been  introduced  and  a  cone  of  length 
50  was  chosen.  Up  to  z  =  10  the  mesh  was  laid  out  in  incre¬ 
ments  of  A<J>  of  15°  with  10  radial  bands  of  mesh  quadrilaterals. 
At  z  =  10,  the  value  of  A<p  was  reduced  to  10°  with  16  radial 
bands  of  mesh  quadrilaterals.  At  z  =  20  the  increments  in  <J> 
were  reduced  to  5°  with  31  radial  bands  of  grid  quadrilaterals. 

In  the  calculated  case  the  solution  was  started  at  small  z 
with  the  a  =  0°  solution  for  a  cone  having  the  9.46°  semi-apex 
angle  and  the  angle  of  attack  was  increased  to  a  =  25°  in 
100  steps.  The  solution  was  continued  downstream  846  steps 
to  a  nominal  distance  of  z  =  50  in  arbitrary  units.  Let  us 
examine  the  flow  field  about  the  cone  at  z  =  50.  If  we  plot 


particle  paths*  in  the  crossflow  plane  as  shown  in  Figure  9(a), 
they  will  not  be  tangent  to  the  cone  since  the  conical  surface 
is  expanding.  If  we  transform  from  Cartesian  coordinates  to 
spherical  polar  coordinates,  this  problem  disappears  since  the 
streamlines  are  conical  as  seen  by  an  observer  at  the  apex. 

Figure  9(b)  shows  this  picture  based  on  the  flow  field  calculated 
at  z  =  50.  It  is  seen  that  the  streamlines  are  altogether 
different  from  the  particle  paths  in  the  crossflow  plane.  It 
appears  that  separation  occurs  near  the  top  of  the  body. 

The  occurrence  of  separation  with  the  use  of  the  Euler 
equations  has  been  observed  many  times.  It  is  usually  attri¬ 
buted  to  artificial  viscosity  due  to  the  numerics.  It  is, 
however,  interesting  to  see  how  the  separation  angle  obtained 
by  the  Euler  equations  with  no  Kutta  conditions  compares  with 
experiment.  The  calculated  separation  points  are  shown  in 
Figure  10  as  a  function  of  axial  length  along  the  cone.  It 
is  seen  that  no  separation  occurs  until  z  >  20,  and  the  separa¬ 
tion  travels  down  the  side  of  the  body.  The  solution  is  not 
valid  since  the  separation  line  is  not  conical  for  the  whole 
length  of  the  cone  as  found  experimentally. 

It  was  decided  to  carry  out  the  calculations  to  larger 
values  of  z  to  see  if  the  separation  would  approach  a  conical 
state  asymptotically.  The  calculations  were  extended  to  z  =  250 
where  <J> s was  approximately  150°,  but  there  was  no  firm  indica¬ 
tion  that  an  asymptotic  solution  was  being  approached.  The 
calculations  were  not  calculated  further  because  of  cost.  It 
is  possible  that  an  asymptotic  separation  position  might  have 
been  achieved.  However,  there  is  no  a  priori  reason  to  believe 
that  the  value  of  4>s^  would  agree  with  experiment.  This  is  to 
be  expected  in  view  of  the  fact  that  the  artificial  viscosity 
in  the  calculation  cannot  be  expected  to  be  equal  to  that  for 
a  turbulent  boundary  layer.  The  Euler  equations  will  not  auto¬ 
matically  yield  the  proper  separated  solution  if  left  to  their 

*Partical  paths  are  solutions  to  the  equations  dx/dy  =  u/v. 
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own  devices.  Instead  we  must  impose  separation  where  we  want  it, 
and  then  apply  the  Euler  equations.  This  procedure  is  essential 
in  obtaining  efficient  starting  solutions  at  high  angles  of 
attack  for  configurations  with  pointed  noses. 

5.3  Solution  with  Kutta  Condition 

The  Euler  equation  applied  to  a  cone  will  yield  a  conical 
solution  if  the  cone  is  not  above  its  shock  detachment  angle. 

If  we  specify  conical  separation  boundary  conditions,  then  the 
solution  should  still  be  conical  since  no  characteristic  length 
has  been  introduced  into  the  problem.  For  an  angle  of  attack  of 
25°,  a  Mach  number  of  3.0,  and  a  cone  angle  of  9.46°  we  find 
that 

ip  =  2.64  Mn  =  1.27 
From  Figure  6  we  find  that 


«t>-  -  132° 

S1 

In  the  calculation,  we  have  used  an  increment  in  A<J>  of  5° 
and  have  used  <|>Sl  =  130°  to  put  separation  on  a  nodal  point 
of  the  mesh.  Let  us  now  address  the  initial  conditions  and 
the  separation  boundary  conditions. 

The  initial  condition  was  taken  to  be  the  solution  for 

calculation  with  no  Kutta  condition  at  z  =  50  where  <J>_  =  162°. 

S1 

The  calculation  was  initialized  at  z  =  0.1  on  this  basis  and 
the  Kutta  condition  was  turned  on  at  z  =  0.2.  In  the  Kutta 
condition  we  have  used  the  conditions  that 

♦a  "  6c 
*c  -  90° 

so  that  Equation  (3)  yields 
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(6) 


u  =  q  sinQ  sin<b 
o  c  Ys 

v„  =  q  sin8  cos* 
o  no  c  rs 

w^  =  q  cos0„ 
o  no  c 

•>  — ► 

Letting  i,  j,  and  k  be  unit  vectors  along  x,  y,  and  z, 
respectively,  and  remembering  that  x,y,z  is  a  left-handed 
Cartesian  system,  we  find  that  the  unit  vector  lying  in  the 
cone  along  the  <)>  =  <j>  line  is 

— ►  -► 

ir  =  i  sinecsin4>s  +  j  sin0ccos<t>s  +  k  cos0c  (7) 

It  is  clear,  therefore,  that  qQ  lies  along  the  ir  direction. 
With  this  version  of  the  Kutta  condition,  the  calculations 
were  carried  out  on  a  grid  with  a  A<J>  increment  of  5°  and  with 
31  rings  in  the  r  direction.  The  calculations  were  carried  to 
z  =  10  without  changing  the  mesh. 

In  Figure  11  the  calculated  flow  field  streamlines  as 
seen  by  an  observer  at  the  cone  apex  are  shown.  The  flow  does 
separate  at  the  $  *  130°  lines,  shown  on  the  figure,  and  no 
secondary  separation  is  apparent.  It  is  of  interest  to  show 
the  degree  of  conicity  of  the  surface  pressure  distribution, 
the  bow  shock,  and  the  flow  field.  Also  the  comparison  between 
the  experimental  and  calculated  pressure  distribution  will  be 
shown.  Some  account  of  the  magnitude  and  sources  of  the 
vorticity  of  the  flow  field  is  then  given. 

The  calculated  body  pressure  coefficients  at  z  =  2 ,  4 , 

6,  8,  and  10  are  shown  in  the  following  table.  The  body 
pressure  distribution  exhibits  a  high  degree  of  conicity. 
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m 


±_ 

z  =  2 

z  =  4 

z  =  6  z  =  8 

z  =  10 

0 

0.6159 

0.6158 

0.6158  0.6158 

0.6159 

30 

0.4997 

0.4995 

0.4993  0.4992 

0.4991 

60 

0.2451 

0.2447 

0.2445  0.2444 

0.2443 

90 

0.0180 

0.0183 

0.0183  0.0183 

0.0183 

120 

-0.1068 

-0.1069 

-0.1069  -0.1070 

-0.1065 

150 

-0.1129 

-0.1116 

-0.1101  -0.1098 

-0.1098 

180 

-0.0678 

-0.0680 

-0.0686  -0.0688 

-0.0689 

The 

bow  waves , 

if  conical. 

will  exhibit  constant 

value  of 

Rg  for  a 

given  <}>  as 

z  increases. 

The  following  values 

are 

given 

R 

R 

s 

s 

z  <(>  =  0 

<j)  =  180 

2  1.584 

5.700 

4  1.584 

5.753 

6  1.584 

5.778 

8  1.584 

5.794 

10  1.584 

5.806 

The  bow  wave  is  closely  conical. 

The  flow  field  on  the  horizontal  plane  of  symmetry  of 
the  body  (<j>  =  90°)  is  now  exhibited. 


o 

o 

z  =  4.0 

z  =  10 

R 

U 

V 

W 

U 

V 

W 

1.0000 

0.  xll33 

0.37369 

0.67022 

0.11133 

0.37371 

0.67022 

1.6387 

0.14217 

0.28260 

0.67055 

0.14216 

0.28256 

0.67056 

2.2794 

0.14062 

0.26327 

0.67224 

0.14064 

0.26324 

0.67223 
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As  a  second  sample  of  the  flow  field,  consider  a  line,  <p  =  130°, 
through  the  separation  point. 


z  =  4.0 

z  =  10 

R 

U 

V 

W 

U 

V 

W 

1.0000 

0.10576 

0.08874 

0.82833 

0.10575 

0.08873 

0.82827 

1.0880 

-0.10190 

0.45520 

0.70642 

-0.10128 

0.45388 

0.70646 

1.1760 

-0.10610 

0.43003 

0.70588 

-0.10572 

0.43033 

0.70571 

3.6401 

0.05600 

0.34271 

0.69840 

0.05589 

0.34266 

0.69850 

Both  samples  of  the  flow  field  show  that  it  is  closely  conical. 

It  is  of  interest  to  see  the  calculated  pressure  distribu¬ 
tion  at  the  cone  surface  since  it  was  essentially  conical.  The 
calculated  distribution  for  the  present  example  is  shown  in 
Figure  12  as  plotted  against  conical  polar  angle.  The  calcu¬ 
lated  pressure  distribution  is  smooth  on  the  windward  surface 
and  on  the  leeward  surface  up  to  the  separation  point.  There¬ 
after  it  is  no  longer  smooth  as  a  result  of  shock  and  expansion 
waves.  The  pressure  coefficient  almost  approaches  Cpmin>  the 
value  for  a  vacuum. 

The  pressure  coefficient  has  been  plotted  versus  lateral 
distance  from  the  vertical  plane  of  symmetry  in  Figure  13. 

The  pressure  coefficient  is  shown  for  the  calculation  without 
the  Kutta  condition  and  with  the  Kutta  condition.  Also 
included  on  the  figure  are  the  measured  pressure  distributions 
from  the  results  of  Landrum,  Reference  22.  It  is  clear  that  the 
data  are  in  close  agreement  with  the  calculation  based  on  the 
Kutta  condition  except  at  <f>  *  180°.  The  use  of  the  Kutta 
condition  changes  the  pressure  distribution  only  on  the  leeward 
side  above  the  separation  point. 

Let  us  now  examine  the  circulation  for  several  contours 
in  the  crossflow  plane.  The  first  contour  is  one  which  enclose 
the  region  between  the  bow  wave  (including  the  bow  wave)  and 
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the  body  contour  as  shown  in  Figure  14  constructed  from  the 
calculations  for  the  case  with  no  Kutta  condition.  Up  to 
z  =  20  there  is  no  separation  yet  a  circulation  is  calculated. 
This  could  be  associated  with  the  starting  solution.  However, 
for  z  >  20  the  circulation  parameters  start  to  increase  more 
rapidly  as  would  be  expected  for  flow  with  separation. 

More  to  the  point  is  the  circulation  as  calculated  for 
the  case  with  the  Kutta  condition.  The  total  circulation  is 
shown  in  Figure  15  by  the  circles,  and  it  increases  almost 
linearly  as  would  be  expected  for  conical  flow.  A  separate 
contour  was  used  enclosing  only  the  bow  shock  to  determine  the 
magnitude  of  the  vorticity  generated  by  it  as  measured  by  the 
circulation  about  the  contour.  It  turns  out  to  be  quite  small, 
about  two  orders  of  magnitude  less  than  that  due  to  separation, 
and  is  difficult  to  calculate  accurately.  The  amount  of 
circulation  below  the  <f>  =  90°  line  was  also  calculated  to  see 
how  much  of  the  total  vorticity  is  on  the  windward  side.  The 
quantity  could  not  be  accurately  determined  but  it  is  about 
1 %  to  2 %  of  the  total  vorticity. 

The  foregoing  calculation  has  demonstrated  the  method  for 
obtaining  solutions  for  the  flow  about  pointed  cones  at  angles 
of  attack  above  the  cone  semi-apex  angle  but  below  the  cone 
shock  detachment  angle.  These  solutions  are  required  for 
starting  Euler  solutions  for  supersonic  configurations  with 
sharp  noses.  It  is  believed  the  separated  flow  solutions  are 
accurate  representations  of  the  flow  field  when  secondary 
separation  does  not  occur  or  when  it  is  not  important.  Further 
calculations  are  needed  to  determine  the  limitations  of  the 
method  with  respect  to  cone  angle,  angle  of  attack,  and  Mach 
number. 


i 
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6.  SAMPLE  CALCULATIONS  FOR  WING-BODY  COMBINATIONS 


6.1  Introductory  Remarks 

It  is  our  desire  to  do  a  complete  pressure  distribution 
calculation  for  a  sample  wing-body  combination  to  illustrate 
the  various  nonlinear  aerodynamic  phenomena  influencing  the 
loading  and  to  compare  with  experiment  where  possible.  For 
this  purpose  we  have  chosen  the  following  configuration.  The 


case  of  M^  =  3.0  and  a  =  25°  will  be  taken.  We  will  compare 
pressure  distributions  and  axial  loading  distributions  on  the 
body  as  measured  and  as  predicted  up  to  the  beginning  of  the 
winged  part.  We  will  compare  the  fin  normal  force  with 
measurements  on  a  similar  body  using  a  fineness  ratio  three 
ogival  nose  rather  than  a  conical  nose.  This  difference  should 
not  be  significant. 

The  solution  for  the  conical  nose  in  the  previous  section 
was  used  to  start  the  solution.  For  z/D  =  3  to  z/D  =  7  a  mesh 
of  36  azimuthal  points  (every  A<J>  =  5°)  was  used  with  30  radial 
points.  From  z/D  =  7.0  to  z/D  =  10  the  number  of  radial  points 
was  increased  to  50.  From  z/D  =  7  to  z/D  =  8 ,  an  integration 
interval  of  Az/D  =  0.004  was  used,  and  this  interval  was  halved 
over  the  winged  section  of  the  combinations. 

On  the  cone  the  separation  parameters  of  the  previous 
section  were  used. 


130° 


90° 


9.46° 


= 


K  - 


On  the  cylinder  the  following  separation  parameters  were  used. 


4>s  =  120°  4>c  =  90°  4>a  =  0° 

At  the  fin  leading  edge  the  velocity  vector  was  taken  tangent 
to  the  leading  edge  in  fulfillment  of  the  Kutta  condition. 

In  the  following  discussions  we  will  first  discuss  the 
pressure  distributions  and  loading  on  the  body  up  to  the  wing 
section  and  then  discuss  the  results  for  the  winged  section. 


6.2  Body  Pressure  Distributions  and  Loadings 


Before  discussing  the  pressure  distributions  let  us  look 
at  the  axial  variation  of  the  normal  loading.  This  can  be 
most  simply  considered  in  terms  of  an  effective  crossflow  drag 
coefficient. 


Cpr  cos<)>  d<J> 


.  ? 

D  sm  a 


(8) 


r  =  local  body  radius 
D  =  body  diameter 

The  appearance  of  D  in  the  denominator  normalized  c^  to  the 

2  ^ 

cylindrical  diameter  and  sin  a  put  Cdc  on  the  basis  of  cross- 
flow  dynamic  pressure.  From  the  Euler  calculation,  crossflow 
pressure  distributions  at  a  number  of  axial  stations  have  been 
determined  and  have  been  integrated  to  obtain  the  effective 
crossflow  drag  coefficient.  The  results  are  shown  in  Figure  16. 

Up  to  z/D  =  3,  the  flow  field  is  conical  so  that  cdc,  which 
used  the  cylinder  diameter  as  the  reference  length,  increases 
linearly  as  shown.  However,  crossing  the  nose-cylinder  juncture 
the  loading  falls  discontinuously .  On  the  upper  meridian 
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(<p  =  180°)  and  the  lower  meridian,  the  flow  expands  around  the 
juncture  by  9.46°.  However,  the  expansion  on  the  upper 
meridian  starts  from  a  much  higher  Mach  number  than  that  on  the 
lower  meridian,  so  that  it  does  not  fall  so  far  in  pressure. 

For  instance,  the  average  upper  surface  pressure  coefficient 
is  -0.10  before  expansion  and  -0.135  after  expansion.  For 
the  lower  surface  the  comparable  pressure  coefficients 
are  0.45  and  0.19.  From  z/D  =  3.0  to  z/D  =  6.0,  there  is 
not  much  change  in  c,jc  ,  its  value  falling  from  about  1.8  to 
about  1.55  in  three  body  diameters. 

Pressure  distributions  are  available  at  a  number  of  body 
stations  from  Reference  22.  In  fact,  the  conical  pressure 
distributions  calculated  for  the  conical  nose  have  already 
been  favorably  compared  with  data  in  Figure  12.  Pressure  dis¬ 
tributions  as  calculated  and  measured  are  compared  for  5  axial 
stations  on  the  cylinder  in  Figure  17.  Examination  of  these 
comparisons  shows  that  the  pressure  coefficients  on  the  impact 
surface  are  predicted  as  closely  as  could  be  expected.  (The 
data  are  for  M  *  2.98  and  the  theory  is  for  M  =  3.0.)  The 
theory  on  the  upper  surface  indicates  a  shock  in  front  of  the 
separation  and  an  expansion  behind  it,  which  is  not  borne  out 
by  the  data.  Switching  to  a  finer  grid  at  z/D  =  7  reduced  this 
shock  jump  from  about  0.1  at  z/D  =  6.5  to  about  0.06  at 
z/D  =  8.  It  thus  is,  in  part,  grid  dependent.  Other  factors 
are  probably  also  responsible.  However,  the  average  pressure 
on  the  upper  surface  is  quite  well  represented  by  the  calculated 
pressure  distribution.  It  follows  that  the  local  loading  is 
closely  predicted  by  the  Euler  code.  The  only  exception  is 
possibly  at  z/D  =  3.8  where  the  measured  loading  is  about  12 
percent  low.  This  is  thought  to  be  a  viscous  effect  associated 
with  the  sharp  corner  at  the  cone-cylinder  juncture. 

It  has  already  been  stated  that  only  primary  separation 
has  been  involved  in  the  present  model.  The  inclusion  of 
secondary  separation  in  the  model  is  possible  but  is  not 
warranted  on  the  basis  of  the  loading  results.  It  is  not  clear 
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that  the  neglect  of  secondary  separation  is  the  reason  for 
the  pressure  fluctuations  on  the  upper  surface  since  these 
seem  to  be  mesh  dependent. 

6.3  Aerodynamics  of  the  Finned  Section 
of  the  Combination 

Several  interesting  aspects  of  the  aerodynamics  of  the 
finned  section  of  the  wing-body  combination  are  shown  by  the 
Euler  calculations.  First,  let  us  examine  the  pressure  dis¬ 
tributions  acting  on  the  combination  in  crossflow  planes 
normal  to  the  axis.  These  are  shown  in  Figure  18  for  fin 
axial  location  starting  at  z/D  =  8.0,  the  start  of  the  fin 
root  chord,  to  z/D  =  10  at  the  fin  trailing  edge.  At  z/D =8.0 
the  upper  surface  pressure  shows  a  crossflow  shock  in  front 
of  separation,  the  magnitude  of  which  is  mesh  dependent  as 
discussed  in  connection  with  Figure  17.  However,  this  shock 
and  the  subsequent  expansion  damp  out  quickly  over  the  finned 
section,  and  the  upper  surface  pressure  becomes  quite  uniform 
by  the  time  the  trailing  edge  is  reached. 

On  the  lower  surface  there  is  a  compression  where  the 
body  fairs  into  the  fin  followed  by  a  later  expansion  slightly 
before  the  fin  tip  is  reached.  There  is  some  roughness  in 
the  calculation  on  the  curved  undersurface  of  the  wing. 

Since  we  have  included  realistic  body  vortex  effects  in 
the  present  calculation,  the  fin  normal-force  coefficient  as 
calculated  should  be  comparable  with  experiment.  A  numerical 
integration  of  the  fin  pressure  coefficients  yields  a  normal- 
force  coefficient  for  the  fin  based  on  fin  planform  area  of 
0.468.  An  experimental  value  of  the  fin  normal-force  coef¬ 
ficient  from  Reference  21  gives  0.457.  This  agreement  is  as 
good  as  can  be  expected  considering  experimental  ^rror  and 
the  fact  that  the  thickness  effects  are  not  included  in  the 
Euler  calculations. 

Some  idea  of  the  character  of  the  vortex  field  and  the 
loading  on  the  body  due  to  the  fin  is  given  by  the  Euler 


calculations.  At  the  beginning  of  the  finned  section,  the 
Kutta  condition  at  the  body  ( 0 s ^  =  120°)  is  no  longer  carried 
along  in  the  calculation.  However,  what  happens  to  the  body 
vortex  in  view  of  the  persistency  of  vorticity  is  of  interest. 
Examination  of  the  crossflow  velocity  components  on  the  body 
yields  an  approximate  idea  of  the  location  of  separation  since 
the  body  is  almost  cylindrical,  the  effects  of  blending  being 
small  near  the  separation  lines.  The  lateral  positions  of  the 
separation  point  were  estimated  to  vary  with  z/D  as  given  below. 


z/D 

xs/D 

8.0 

0.44 

8.25 

0.43 

8.50 

0.44 

8.75 

0.48 

9.00 

0.48 

9.25 

0.47 

9.50 

0.46 

9.75 

0.44 

10.0 

0.43 

Since  a  value  of  xs/D  of  0.444  corresponds  to  a  separation  point 
at  4>s j  =  120°,  it  is  seen  that  the  separation  line  varies  above 
and  below  the  <j>Sl  =  120°  line.  Within  the  accuracy  of  the 
determination  of  xs/D,  it  is  well  approximated  by  a  line  par¬ 
allel  to  the  body  centerline  over  the  winged  part  of  the 
configuration.  However,  the  persistency  of  circulation  has  per¬ 
petuated  separation  without  imposing  the  Kutta  condition. 

Since  the  body  vortices  persist  over  the  fins,  what  can 
we  say  about  the  "lift  carry-over"  from  the  fin  to  the  body? 
Figure  19  shows  the  effective  crossflow  drag  coefficient  curve 
of  Figure  16  extended  over  the  fin.  It  is  noted  that  a 
linear  rate  of  increase  in  c^  starts  at  the  leading  edge  of 
the  fin  root  chord.  This  rate  of  increase  is  composed  of 
components  associated  with  the  average  pressure  Cp  on  the 


upper  surface  and  the  average  pressure  Cp^  on  the  lower  sur¬ 
face.  It  turns  out  that  CpUp  does  not  change  signif icantly 
going  from  the  cylinder  to  the  finned  surface.  (The  average 
pressure  is  almost  that  of  a  vacuum.)  Accordingly,  the  body 
lift  on  the  upper  surface  is  not  changed  by  fin-body  inter¬ 
ference,  but  is  the  same  as  that  in  front  of  the  fin  with  the 
body  vortices  present.  The  increase  in  c<ic  over  the  finned 
section  is  due  mainly  to  an  increase  in  Cp^.  This  effect  is 
associated  with  the  damming  effect  of  the  fin  on  the  flow  on 
the  impact  side  of  the  configuration. 


7.  EFFECTS  OF  ANGLE  OF  ATTACK  AND  MACH  NUMBER  ON  LOADING 
OF  EXAMPLE  WING- BODY  COMBINATION 

In  order  to  demonstrate  the  use  of  the  Euler  code  for 
determining  pressure  distributions,  it  was  decided  to  do  some 
parametric  calculations  for  a  sample  wing-body  calculation 
varying  angle  of  attack  and  Mach  number.  The  basic  configura¬ 
tion  chosen  for  these  calculations  has  the  following  dimensions: 


The  configuration  has  a  cruciform  wing  and  an  ogival  nose.  At 
each  axial  station  of  the  computation  the  cross  section  of  the 
configuration  is  transformed  analytically  into  a  circle  of 
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unit  radius  assuming  the  fins  have  no  thickness.  The  boundary 
conditions  are  applied  on  a  circle  of  radius  1  +  e  where 
e  =  0.025  at  the  leading  edge  of  the  wing-body  juncture.  The 
value  of  e  varies  from  0  at  z  =  4.16  to  0.025  at  z  =  5.16. 
Since  the  boundary  conditions  are  applied  slightly  off  the 
unit  circle,  the  wing  is  faired  smoothly  into  the  body  and 
forms  a  blended  wing-body  combination.  The  wing  has  thickness 
as  a  result. 

The  separation  line  position  on  the  body  was  not  varied 
during  the  calculations  since  a  strictly  parametric  study  was 
desired.  The  separation  line  was  taken  at  <{>Sl  =  120°  from 
z  =  1.5  to  z  =  5.16.  The  separation  angles,  defined  in 
Figure  2,  were  taken  to  be 

<J>c  =  30°;  <j)a  =  10° 

The  streamline  at  separation  was  thus  in  a  plane  30°  from  the 
tangent  plane  to  the  body  and  10°  from  the  axial  direction. 
This  prevented  it  from  entering  the  body  on  the  ogival  nose. 
These  conditions  were  kept  constant  for  the  following  set  of 
calculations . 


M  =  3  a  =  15°,  20°,  25° 

OO  99 

a  =  20°  =  3,  4,  5 

A  variable  grid  was  used  in  calculating  the  separated 
flow  up  to  the  leading  edge  of  the  wing-body  juncture.  How¬ 
ever,  over  the  winged  part  of  the  configuration,  a  constant 
mesh  was  used  with  A^>  =  3°  and  51  radial  rings  of  quadri¬ 
laterals.  The  integration  along  the  wing  was  carried  out  in 
1000  steps  of  equal  axial  distance  (0.002). 

It  is  of  interest  to  examine  first  the  effect  of  Mach 

number  on  the  fin  pressure  distributions  at  the  wing  trailing 

edge  as  shown  in  Figure  20(a).  On  the  impact  surface  the 

2 

pressure  is  not  far  from  the  Newtonian  value  of  2  sin  a  and 
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does  not  vary  significantly  with  Mach  number.  On  the  upper 
surface  the  pressure  coefficient  is  nearly  uniform  and  is 
close  to  that  for  zero  pressure,  a  value  dependent  on  Mach 
number.  By  using  the  above  approximations,  the  fin  normal 
force  could  be  estimated. 

At  Mm  =  5,  the  body  upper  surface  pressure  coefficient  at 

the  trailing  edge  varies  from  -0.042  to  0.0498  to  be  compared 

with  a  value  of  Cp  .  of  -0.0574.  At  a  Mach  number  of  3.0  the 

“mm 

body  upper  surface  pressure  coefficient  varies  from  about  -0.10 

to  -0.12  versus  a  value  of  Cp  -  of  -0.159. 

Miun 

Comparison  of  the  fin  pressure  distributions  at  the  fin 
trailing  edge  at  =  3  for  angles  of  attack  15°  and  25°  is 
made  in  Figure  20(b).  The  pressures  on  the  upper  surfaces 
are  close  to  Cp  .  ,  the  a  =  25°  pressures  being  closer  than 
the  a  =  15°  ones.  On  the  impact  surface  the  pressure  shows 
considerable  difference  between  a  =  15°  and  a  =  25°  as  would 
be  expected  on  the  basis  of  Newtonian  impact  theory. 

A  point  of  interest  in  connection  with  the  foregoing  cal¬ 
culations  is  how  the  vorticity  due  to  the  body  compares  in 
magnitude  with  that  due  to  the  fin  leading  edge.  The  circula¬ 
tion  about  a  contour  similar  to  that  in  Figure  14  is  useful 
in  illustrating  these  effects.  The  value  of  r/2-nVma  has  been 
plotted  versus  axial  distance  in  Figure  21(a)  for  a  =  15° 
and  25°  at  M^  =  3.0.  The  amount  of  vorticity  generated  before 
the  separation  point  on  the  body  is  negligible,  but  on  the 
cylindrical  section  it  grows  approximately  linearly  until  the 
winged  section  is  reached.  Over  the  winged  part  of  the  missile, 
the  rate  of  vorticity  generation  is  smaller  at  a  =  15°  than  at 
a  =  25°  as  would  be  expected. 

A  comparison  of  the  circulation  distributions  for  Mm  =  3 

and  5  for  a  =  20°  is  given  in  Figure  21(b).  Over  the  body  up 

to  the  winged  section,  the  circulation  results  show  a  'ightly 
lower  value  for  M  =  5.0  than  for  M  =  3.0,  but  the  differences 

are  not  large.  Over  the  winged  part  of  the  body  it  would  be 
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5  than 


expected  that  less  vorticity  would  be  generated  at  = 

at  M  =3  since  the  fin  normal  force  is  somewhat  less  (see 
00 

Fig.  20(a)). 

8.  COMPARISON  OF  BODY  VORTEX  EFFECTS  PREDICTED  BY 
EULER  EQUATIONS  AND  A  PANEL  METHOD 

8.1  Introductory  Remarks 

It  is  possible  to  obtain  the  fin  and  body  pressure  dis¬ 
tributions  with  or  without  body  vortices  present  using  either 
the  Euler  equations  or  a  supersonic  panel  method  such  as  that 
given  in  Reference  2.  The  panel  method  is  basically  a  linear 
method  for  the  velocity  flow  field,  even  though  it  uses  the 
nonlinear  Bernoulli  equation  for  the  pressure  field.  We  have 
chosen  as  the  example  the  following  wing-body  combination  at 
a  =  10°  and  M^  =  3.0.  The  comparison  between  the  Euler  method 
and  the  panel  method  will  thus  be  in  the  possible  range  of 
applicability  of  the  panel  method. 

The  wing-body  combination  has  a  fineness  ratio  3  ogival 
nose  followed  by  a  cylindrical  section  6  body  diameters  long 
up  to  the  winged  section.  The  -  1  delta  fin  has  a  root 
chord  of  2  diameters  and  an  exposed  semispan  of  0.5  diameter. 


For  the  purpose  of  comparing  body  loadings  the  same  body 
vortex  model  is  used  for  both  the  Euler  and  panel  methods.  It 


was  assumed  that  separation  occurred  from  z/D  =  3  to  z/D  =  9 
along  the  4>Sl  =  90°  line.  The  angular  parameters  at  separation 
were 

*a  =  10° 

4C  -  30° 

for  the  entire  cylindrical  length. 

In  the  panel  method  the  vortex  downwash  distribution  across 
the  wing  panels  is  determined  from  external  and  image  line 
vortices  using  the  Biot-Savart  law,  thereby  neglecting  com¬ 
pressibility.  First  we  will  compare  the  vortex  induced  down- 
wash  as  calculated  by  the  Euler  equation  and  then  the  Biot-Savart 
law.  Next  we  will  determine  the  fin  loading  distributions  due 
to  the  Euler  downwash  distribution  using  both  the  Euler  and 
panel  methods.  It  is  noted  that  the  loading  distribution  cal¬ 
culated  for  the  Euler  equations  is  for  a  fin  with  thickness  and 
that  using  the  panel  method  is  for  a  fin  without  thickness. 

8.2  Comparison  of  Downwash  Distributions 

We  will  compare  the  downwash  distributions  resulting  from 
the  vortical  flow  field  at  x/D  =  9  calculated  by  the  panel 
code  and  by  the  Euler  code  with  the  above  separation  model.  The 
downwash  due  to  the  body  vortices  in  the  plane  of  the  wing 
calculated  by  the  Euler  equations  was  determined  by  subtracting 
the  downwash  calculated  with  no  Kutta  condition  from  that  with 
the  Kutta  condition.  This  procedure  eliminates  the  downwash  that 
occurs  without  body  vortices.  The  body  vortex  downwash  in  the 
plane  of  the  wing  as  used  in  the  panel  method  is  determined  by 
considering  the  body  vortex  strengths  and  positions  in  the  cross- 
flow  plane,  putting  the  appropriate  images  inside  the  body,  and 
calculating  the  induced  velocities  by  the  Biot-Savart  law.  In 
the  present  calculation  there  are  60  mesh  quadrilaterals  in  the 
<Ji  direction  and  49  in  the  radial  direction.  The  vorticity  in 
each  quadrilateral  was  determined  from  the  Euler  solution  by 
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calculating  the  circulation  around  its  boundary,  and  a  point 
vortex  was  placed  at  its  centroid.  The  downwash  on  the  hori¬ 
zontal  plane  was  then  determined  by  considering  all  external 
and  image  vortices  on  each  side  of  the  combination,  calculating 
their  induced  downwash,  and  summing  over  all  the  vortices. 

The  calculated  results  are  shown  in  Figure  22  for  both 
the  Euler  method  and  the  Biot-Savart  method.  The  induced 
flow  is  downward  over  the  entire  fin  span.  The  two  results 
are  fairly  close  although  the  Euler  solution  yields  on  the 
average  more  downwash  than  the  Biot-Savart  method.  Since  the 
Biot-Savart  method  neglects  compressibility,  some  difference 
is  to  be  expected. 

8.3  Comparison  of  Loading  Distributions 

In  order  to  make  a  meaningful  comparison  of  the  loading 
distributions  calculated  by  the  Euler  equations  and  the  panel 
method,  it  was  decided  to  use  the  same  downwash  distribution 
at  z/D  =9.0  for  both.  The  loading  is  defined  to  be  the 
lower  surface  pressure  coefficient  minus  the  upper  surface 
pressure  coefficient  at  the  same  position  on  the  planform.  By 
comparing  loadings  we  minimize  any  effects  due  to  thickness 
accounted  for  in  the  Euler  approach  but  not  in  the  panel  method 
since  to  the  second  order  in  a,  thickness  has  no  effect  on  the 
normal-force  coefficient  for  supersonic  airfoils. 

The  comparisons  are  presented  in  Figure  23.  In  Figure  23(a) 
the  Euler  pressure  distributions  in  the  upper  and  lower  sur¬ 
faces  are  presented  with  and  without  the  Kutta  condition.  It  is 
noted  that  the  effects  of  body  vortices  are  larger  on  the 
upper  surface.  However,  near  x/D  =  0.5  on  the  lower  surface, 
the  position  of  the  shock  is  changed  in  a  way  which  is  probably 
a  function  of  the  algorithm  sensitivity.  In  Figure  23(b)  the 
loading  as  determined  from  the  pressure  distributions  in 
Figure  23(a)  is  shown.  The  loading  at  the  wing-body  juncture 
from  the  Euler  solution  may  be  spurious  because  of  the 
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change  in  the  shock  location  noted  above.  Also  in  Figure  23(b) 
the  loading  as  computed  by  the  panel  method  is  compared  with 
that  due  to  the  Euler  code.  There  are  differences,  but  the 
average  load  is  about  the  same. 

The  panel  method  utilized  as  a  downwash  boundary  condition 
on  the  fins  of  zero  thickness  the  difference  in  downwash  between 
the  Euler  code  with  and  without  the  Kutta  condition.  As  a 
matter  of  interest  the  vorticity  axial  distribution  for  the 
body  alone  up  to  z/D  =  9.0  with  and  without  the  Kutta  condi¬ 
tion  is  shown  in  Figure  24(a).  The  nonzero  values  of  r/2-iTV^a 
for  the  body  without  the  Kutta  condition  are  not  due  to 
body  flow  separation  since  no  separation  due  to  numerical 
viscosity  occurred  before  z/D  =  9.0.  The  downwash  in 
of  the  fin  due  to  the  vorticity  distribution  without  the  Kutta 
condition  causes  very  small  downwash  angles  relative  to  those 
with  the  Kutta  condition.  The  panel  calculations  were  made 
with  10  panels  chordwise  and  10  rows  of  panels  spanwise. 

There  are  a  number  of  reasons  why  the  panel  method  is 
only  approximate.  First  the  loading  calculations  were  made 
at  a  =  0  using  the  fin  boundary  conditions  described  above. 

We  could  have  calculated  the  loading  at  a  =  10°  with  and  with¬ 
out  body  vortices  and  subtracted  the  differences.  The  results 
would  be  different  from  those  in  Figure  23  if  the  full  nonlinear 
Bernoulli  equations  were  used.  However,  it  should  be  noted 
that  the  results  of  Figure  23  were  calculated  with  the  non¬ 
linear  Bernoulli  equation  including  the  effects  of  a  strong 
sidewash  due  to  the  vorticity  distribution.  Without  including 
this  sidewash  effect  the  agreement  is  not  nearly  so  good. 

In  the  second  place  there  are  thickness  effects  included 
in  the  Euler  results  but  neglected  in  the  panel  results. 

However,  by  using  loadings  rather  than  pressure  coefficients 
this  effect  is  minimized.  Thirdly,  the  panel  method  assumes 
that  the  downwash  distribution  at  z/D  =  9  is  the  same  for  all 
z/D  values  over  the  wing.  Both  of  these  latter  effects  could 
be  treated  by  the  panel  method. 
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This  example,  taken  at  the  edge  of  the  range  of  appli¬ 
cability  of  panel  theory,  shows  that  the  method  for  handling 
vortex  effects  in  that  theory  is  basically  sound  for  gross 
forces.  However,  several  nonlinear  effects  which  reduce  the 
accuracy  of  panel  theory  have  been  noted. 


9.  CONCLUDING  REMARKS 

The  primary  purpose  of  this  report  has  been  to  demonstrate 
the  applicability  of  the  Euler  equations  to  the  calculation  of 
flow  about  supersonic  wing-body  combinations  at  high  angles  of 
attack  for  cases  where  significant  effects  of  compressibility 
and  flow  separation  are  both  present.  The  treatment  of 
separated  flows  by  the  Euler  equations  requires  a  knowledge  of 
the  positions  of  the  separation  lines.  Special  boundary  con¬ 
ditions  are  used  on  these  lines  to  make  the  Euler  solution 
unique.  It  was  found  that  methods  such  as  those  of  Stratford 
for  predicting  the  separation  line  location  are  not  valid  at 
high  crossflow  Mach  numbers  since  separation  does  not 
necessarily  occur  in  a  rising  pressure  field.  It  was  thus 
necessary  to  rely  principally  on  empirical  correlation  for 
separation  line  positions. 

Since  most  high-speed  missiles  (and  airplanes)  have  sharp 
noses  which  can  be  fitted  by  a  cone  and  since  the  flow 
separates  over  cones  at  angles  of  attack  in  excess  of  the 
cone  half  angle,  it  is  necessary  to  have  solutions  for  cones 
with  separation  to  start  the  Euler  solutions  at  high  angles 
of  attack.  Accordingly,  the  separation  data  on  cones,  which 
exhibit  conical  separation,  were  correlated  so  that  the  primary 
separation  line  locations  could  be  estimated.  At  the  same 
time  primary  separation  angles  for  two-dimensional  circular 
cylinders  and  circular  cylinders  at  angles  of  attack  were 
correlated  for  application  to  the  cylindrical  section  of  wing- 
body  combinations. 
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The  Euler  equations  were  applied  to  a  cone  at  an  angle 
of  attack  well  above  its  semiapex  angle  for  which  separation 
is  to  be  expected.  Without  applying  the  separation  boundary 
condition,  the  solution  finally  exhibits  nonconical  separation 
as  a  result  of  artificial  viscosity.  Such  a  solution  does  not 
agree  with  experiment.  By  imposing  a  conical  separation 
boundary  condition,  the  Euler  equations  developed  a  conical 
separation  flow  for  which  the  pressure  distribution  was  in 
good  agreement  with  experiment.  This  solution  demonstrated 
how  to  develop  a  proper  starting  solution  for  a  wing-body 
combination  when  the  nose  is  pointed  and  separation  starts 
at  the  nose.  Cone  tables  for  angles  of  attack  greater  than 
the  cone  semi-apex  angle  can  be  developed  by  this  technique. 

Having  demonstrated  how  to  start  the  solution  for  bodies 
or  wing-bodies  when  separation  starts  at  the  nose,  we  carried 
out  a  series  of  calculations  to  demonstrate  the  feasibility 
of  using  the  Euler  equations  to  determine  the  flow  field  of  wing- 
body  combinations  with  both  body  vortices  and  leading-edge 
vortices  present.  A  range  of  Mach  numbers  up  to  5  and  a  range 
of  angles  of  attack  up  to  25°  were  covered  with  results  which 
were  in  good  agreement  with  available  data.  No  attempt  to 
determine  the  limits  on  the  method  with  regard  to  angle  of 
attack  was  carried  out.  However,  it  is  felt  that  the  method 
in  its  present  form  should  work  up  to  the  angle  of  shock 
detachment  for  the  nose. 

Since  the  Euler  equations  treat  the  combined  effects  of 
compressibility  and  vortices  on  a  true  nonlinear  basis,  it  was 
decided  to  apply  them  to  a  case  for  comparison  with  a  supersonic 
panel  method  (DEM0N2)  which  linearily  superposes  these  effects. 

A  case  at  M  =  3.0  and  a  =  10°  was  chosen  which  is  near  the 

OD 

limits  of  applicability  of  linear  panel  methods.  It  was  found 
that  the  Euler  equations  predicted  somewhat  more  downwash  on 
the  horizontal  fins  than  did  the  Biot-Savart  law  used  in  the 
panel  method.  For  the  same  downwash  distribution  on  the  fins, 
both  methods  predicted  about  the  same  normal-force  increment 


due  to  the  body  vortices,  but  the  shapes  of  the  loading  dis¬ 
tributions  differed.  Further  comparisons  are  needed  between 
panel  and  Euler  methods  to  determine  the  accuracy  and  short¬ 
comings  of  panel  methods  for  estimating  the  effects  of  vortices 
on  pressure  distributions. 

While  the  present  work  has  demonstrated  the  applicability 
of  the  Euler  equations  to  supersonic  wing-body  combinations 
with  significant  compressibility  and  vorticity  effects,  it  has 
not  addressed  certain  additional  problems  which  arise  for 
wing-body- tai 1  combinations.  One  of  these  problems  is  the 
rolling-up  of  the  fin  wakes  between  the  wing  and  tail,  and 
the  subsequent  effect  on  the  tail  panels.  Another  effect, 
known  to  be  important  at  high  angles  in  certain  cases,  is  the 
vortex  wake  developed  by  the  afterbody  between  the  wing  and 
tail  particularly  as  influenced  by  the  presence  of  the  wing 
vortices.  These  questions  are  susceptible  to  solution  by 
Euler  methods,  but  the  technology  remains  to  be  demonstrated. 
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APPENDIX  A 


METHOD  FOR  SCALING  SEPARATION  LINE  LOCATION 
FROM  OIL- FLOW  PHOTOGRAPHS 

Oil  flow  pictures  are  presented  in  Reference  20  for  a  cone- 
cylinder  at  Mach  numbers  of  1.6,  2.3,  2.96,  and  4.63  at  angles 
of  attack  of  12°,  20°,  28°,  36°,  44°,  52°,  and  60°  for  most  of 
the  Mach  numbers.  Photographs  of  the  oil  flow  patterns  were 
taken  from  the  side  and  above  the  body  so  that  the  near  side 
edge  of  the  body  (<f>  =  90°)  appears  closer  to  the  bottom  of  the 
body  than  to  the  top  of  the  body.  The  approximate  side  view, 
as  seen  by  the  camera,  and  an  end  projection  are  shown  below. 


The  0  =  90°  line  measured  from  the  bottom  of  the  body  is  shown 
by  an  index  mark  on  each  photograph.  It  is  possible  to  scale  the 
distances  a,  b,  and  zg  from  the  photographs  with  a  divider.  The 
values  of  a  and  b  depend  for  their  accuracy  on  the  edge  contrast 
in  the  photograph. 

From  these  scaled  values  in  arbitrary  units  of  length,  it 
is  possible  to  determine  4>Sl,  the  primary  separation  location, 
given  by 


<fc  =  90°  +  6  + 
r  s  ^  s 


(A-  1) 
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With  the  body  radius  given  as 


r 


a  +  b 
2 


(A-2) 


we  find  that 


s  in  6 


b  -  a 
S"  +  a 


(A-2) 


z 

sin  0  =  — p-  -  sin  6  (A-3) 

s  r 

This  gives  all  the  information  necessary  to  get  <J>S^  from 
Equation  (A- 1) . 

The  photographs  of  the  oil  flow  patterns  have  been  read  for 
flow  distances  along  the  body:  x/l  =  0.325,  0.475,  0.675,  and 
0.925.  The  first  station  is  on  the  cone  and  the  last  three 
are  on  the  cylindrical  part  of  the  body.  The  scaled  values  of 
a,  b,  and  zg  are  given  in  Table  2  together  with  the  position  of 
the  line  of  primary  separation.  In  one  instance  the  photograph 
was  not  sufficiently  distinct  to  locate  the  separation  line. 

The  units  of  a,  b,  and  zg  are  arbitrary  since  the  angles  depend 
only  on  length  ratios. 
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LIST  OF  SYMBOLS 


body  radius 

aspect  ratio  of  wing  formed  by  joining  two  opposing 
fins  together  along  their  root  chords 

crossflow  drag  coefficient,  (Cp^  -  CpUp)/<ln 

pressure  coefficient,  (p  -  p^j/q 

average  pressure  coefficient  on  top  of  body 

average  pressure  coefficient  on  bottom  of  body 

pressure  coefficient  for  p  =  0 


body  diameter  of  cylindrical  section 

2y/(y-1) 

length  of  cone-cylinder  body 
M  sina,  normal  Mach  number 

oo  7 

free-stream  Mach  number 
local  static  pressure 
free-stream  static  pressure 

p/% 

q  sin  a 
‘00 

free-stream  dynamic  pressure 

velocity  of  fluid  at  a  separation  line;  see  Figure  2 
radial  coordinate  in  transformed  plane;  see  Figure  1 
value  of  r  at  body  surface 
value  of  r  at  bow  wave 

cylinder  Reynolds  number  based  on  diameter 


R 


r/rb 

value  of  R  at  bow  shuck 


s  semispan  of  fin  measured  from  body  axis 

u,v,w  components  of  local  velocity  along  x,  y,  and  z, 
respectively 

uo,vo’wo  u’  v’  anc*  'J  ('jmPorients  °f  R0 

x,y,z  Cartesian  coordinates  of  physical  space;  see  Figure  1 

xg  x  coordinate  of  separation  line  on  body 

z,r,<(>  cylindrical  coordinates  of  transformed  space;  see 
Figure  1 

Z  rei(*  ‘  7r/2) 

a  angle  of  attack 

y  ratio  of  specific  heats 

T  circulation  about  a  contour  in  a  half-plane  composed 

of  the  body  surface  and  a  contour  containing  the  bow 
wave;  see  Figure  14 

r  .  ,  circulation  about  a  contour  containing  the  shock  in 

snocK  the  right  half-plane 

0C  semiapex  angle  of  cone 

X  fin  taper  ratio 

(j)  polar  angle  in  cylindrical  coordinates;  see  Figure  1 

<j>Sl  polar  angle  of  primary  separation  on  cone  or 

cylinder 

(|>S2  polar  angle  of  secondary  separation 

4 >a  angle  between  qQ  and  wQ;  see  Figure  2 

<J>C  angle  between  plane  of  separation  and  tangent  to 

cylinder  at  separation;  see  Figure  2 
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TABLE  2 .  CONTINUED 


TABLE  2 .  CONCLUDED 


Transformed  Space 


transformation  between  the  physical  space 
:  transformed  curvilinear  space. 


IfTsK 


Figure  2.-  Velocity  components  and  angles  of  the  vortex  sheet 
defining  the  parameters  for  the  Kutta  condition. 
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Figure  4.-  Effect  of  Mach  number  on  position  of  turbulent 
separation  on  circular  cylinder  in  crossflow. 


symbols 


Figure  6.-  Correlation  of  cone  primary  separation  location 


Primary  Separation  Line  Angular  Position, 


Cone  pressure  coefficient 


Figure  12.-  Calculated  surface  pressure  distribution 
on  cone  with  Kutta  condition 
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—  Circulation  contour 
(encloses  bow  wave) 
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/  00 

/  a  =  25° 
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9  .  =  9.46° 
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O  Calculated  from  circulation 
contour 

z,  arbitrary  units 


Figure  14.-  Growth  of  circulation  parameter  along 
cone  with  no  Kutta  condition  imposed. 


Circulation  parameter,  rr~ —  ,  arbitrary  units 


Figure  15.-  Growth  of  circulation  parameters 
along  cone  with  Kutta  condition 
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Pressure  coefficient 


-  w 


jjjjj 


(a)  z/b  =  3.8 


Figure  17,-  Comparison  of  pressure  distribution 
as  calculated  by  the  Euler  code  and  as  measured 
for  example  wing  body  combination. 
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Bod' 
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(a)  z/D  =8.0 

Figure  18.-  Pressure  distributions  in  crossflow  planes 
of  winged  part  of  example  wing-body  combinations. 
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Upper  surface 
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Upper  surface 


(c)  z/D  =  9.0 


Figure  18.-  Continued. 


Figure  18.-  Continued. 


Fm  pressure  coefficient 


Fin  pressure  coefficient 


0-5  0.6  0.7  0.8  0.9  1.0 

Lateral  distance,  x/s 

(b)  Effect  of  angle  of  attack 

Figure  20.-  Pressure  distributions  at  trailing  edge  of  fin 
mounted  on  circular  body  as  calculated  by  an  Euler 

code. 
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(b)  Effect  of  Mach  number 
Figure  21.-  Concluded. 
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SJ 


With  body  vortices 
——  Without  body  vortices 


-0.05 


+0.05 


+0.10 


(a)  Pressure  distributions 
from  Euler  solutions 
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SI  -  1.0 


at  trailing  edge  of  fin. 


Figure  23.-  Pressure  distribution  on  wing-body  combination  at 
fin  trailing  edge  and  fin  vortex  loading  for  a  a  ■  10° 
and  H.  -  3.0. 
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Figure  24.-  Concluded 
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